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ABSTRACT 


Arrival traffic scenarios with 808 flights from 173 airports to Houston George Bush 
International airport are simulated to determine if Call For Release flights can receive a benefit 
in terms of less delay over other flights by scheduling prior to gate pushback (look-ahead in 
time) as opposed to at gate pushback. Call for Release flights are departures that require 
approval from Air Route Traffic Control Center prior to release. Realism is brought to the study 
by including gate departure delay and taxi-out delay uncertainties for the 77 major U. S. airports. 
Gate departure delay uncertainty is assumed to increase as a function of look-ahead time. 
Results show that Call For Release flights from an airport within the freeze horizon (a region 
surrounding the arrival airport) can get an advantage over other flights to a capacity constrained 
airport by scheduling prior to gate pushback, provided the wheels-off time uncertainty with 
respect to schedule is controlled to a small value, such as within a three-minute window. 
Another finding of the study is that system delay, measured as the sum of arrival delays, is 
smaller when flights are scheduled in the order of arrival compared to in the order of departure. 
Because flights from airports within the freeze horizon are scheduled in the order of departure, 
an increase in the number of internal airports with a larger freeze horizon increases system 
delay. Delay in the given scenario was found to increase by 126% (from 13.8 hours to 31.2 
hours) as freeze horizon was increased from 30-minutes to 2-hours in the baseline scenario. 


I. INTRODUCTION 


Automated tactical scheduling of departure traffic for joining a constrained enroute or arrival 
traffic flow is one of the steps needed for integrating arrival, departure, surface and enroute 
operations. Tactical scheduling of departures is accomplished by a Call For Release (CFR) 
procedure, also known as the Approval Request (APREQ) procedure, which requires the Air 
Traffic Control Tower (ATCT) to seek permission from the Air Route Traffic Control Center 
(ARTCC) prior to departure release. Reference 1 describes the Precision Departure Release 
Capability (PDRC) operational concept and system design for automating the CFR process 
using the “outbound” and “inbound” scheduling functions of the Traffic Management Advisor 
(TMA) developed by the National Aeronautics and Space Administration (NASA) and the 
Federal Aviation Administration (FAA). FAA has since renamed TMA as Time-Based Flow 
Management (TBFM). 

PDRC uses TBFM’s Enroute Departure Capability (EDC) function for outbound scheduling 
of departure flights from airports included in the CFR within the ARTCC to meter points within 
the ARTCC. It uses TBFM’s arrival metering called “internal” departure function for scheduling 
departures into inbound arrival streams of TBFM-metered airports when they are departing from 
an airport within the arrival freeze horizon. In addition to these functions, a Coupled Scheduling 


function for connecting EDC schedules generated in one ARTCC with arrival metering in an 
adjacent ARTCC has been available in TBFM versions since 2011. This function binds the 
outbound and inbound tactical departure scheduling. 

This report is focused on inbound scheduling. It augments the study in Ref. 2 that examined 
the benefit of scheduling CFR flights prior to gate pushback compared to upon gate pushback. 
The earlier study (Ref. 2) assumed gate departure time and taxi-out time uncertainties for 77 
major U. S. airports that were characterized based on real-world data. The study scenario 
consisted of 27 hours of traffic from January 2, 2011 to January 3, 2011 with 37,346 flights from 
1,552 airports to 1,357 airports. Arrival capacities at the 77 major U. S. airports were obtained 
from the FAA’s Aviation System Performance Metrics (ASPM) database. Other airports were 
assumed to be unconstrained so that the traffic bound for these airports would not be delayed 
by the arrival scheduler. All flights to all airports were assumed to be CFR flights. Delayed 
departures resulting from gate departure time and taxi-out time uncertainties were allowed to 
increase speed by up to 2.5% to arrive at the slot reserved at gate pushback or prior to gate 
pushback. Statistics of ground and airborne delays were obtained by running 100,000 Monte 
Carlo simulations with a First-Come First-Served (FCFS) arrival scheduler for each of the five 
look-ahead times- from at pushback (zero look-ahead time) to 60 minutes prior to pushback in 
15-minute intervals. The main finding of the study in Ref. 2 is that scheduling should be done at 
gate pushback because a slot assigned prior to gate pushback is usually lost due to the 
combination of gate departure time and taxi time uncertainties. The cost of delay (weighted 
delay), which is measured as the sum of two-times airborne delay and one-time ground delay, 
increases for longer look-ahead times. Without uncertainty, there is a slight advantage for 
scheduling early as delay otherwise taken in the air can be transferred to the ground; transfer of 
delay to the ground decreases the overall weighted delay. 

This study examines 808 flights originating from 173 airports destined for Houston George 
Bush Intercontinental Airport (IAH) within the 27-hour interval from January 2, 2011 to January 
3, 2011, the same period considered in Ref. 2. The 808 flights to IAH are a subset of the flights 
used in the earlier study. The first objective of using this dataset is to see if there is a marked 
advantage for CFR flights when mixed in with non-CFR ones. The second objective is to 
determine total system delay in several scenarios with CFR and non-CFR flights. As in the 
previous study, this study also employs the concept of look-ahead time and introduces the 
concept of freeze horizon for internal departures. Using an updated FCFS scheduler, three 
overall categories of simulations were run: 1) all flights are CFR flights, 2) all flights are non- 
CFR flights and 3) a mix of CFR and non-CFR flights. 

Analysis of the results of these scenarios show that CFR flights have an advantage of early 
scheduling over non-CFR flights when they originate from airports within the freeze horizon, 
provided the wheels-off time uncertainty of the CFR flights can be controlled to a small value 
like within a three-minute window, and the freeze horizon is smaller than the sum of the look- 
ahead time and the flight time to the destination airport. If the gate departure time uncertainty of 
CFR flights is high, the benefit of early scheduling is lost. Reallocation of new slots for missed 
slots increases delay, especially airborne delay. It is best to schedule at gate pushback in such 
situations. System delay results show that scheduling based on estimated time of arrival in the 
order of arrivals leads to less delay compared to in the order of departures. The implication of 
ordering by departures is that an arrival slot that could have been given to a later departure that 
arrives earlier is not available for assignment because an earlier departure had been assigned 
that slot since it was scheduled earlier. This would result in a delay for the later departure. A 
longer freeze horizon causes more flights to be scheduled in the order of departure, and 
therefore it causes the system delay to increase. However, the greatest system delay was 
observed at the longest look-ahead time in the all CFR flights scenario due to increasing gate 
departure time uncertainty with increasing look-ahead time. 

The rest of the technical memo is organized as follows. Section II describes the IAH arrival 
scenario and defines freeze horizon, look-ahead time, CFR and non-CFR flights, external and 
internal airports and the scheduling process. Section Ill provides descriptions of the gate 
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departure delay and taxi-out delay distribution modeling using ASPM Data, Airspace Concept 
Evaluation System (ACES) and the FCFS scheduler. The block diagram in this section outlines 
the procedure for preparing the input for the scheduler using ACES simulated flight time and for 
scheduling and re-scheduling CFR flights, and non-CFR internal departures and external 
departures. Statistics of ground, airborne and weighted delay for the different traffic scenarios 
are discussed in Section IV. System delay results for these scenarios are also discussed in this 
section. Finally, the main findings are stated in Section V, the conclusions section. 


ll. STUDY SCENARIO AND SCHEDULING PROCEDURE 


The following section describes the traffic scenario used in this study and the manner in 
which scheduling is handled for internal, external, CFR and non-CFR flights. 


A. Traffic Scenario 


The scenario described in this section is 
tailored for scheduling departures to join an @ DEN Freeze horizon 
arrival stream at a capacity constrained arrival *s 
airport. The concept of freeze horizon from ge 
TBFM literature (see Ref. 3) is employed to [225] 
distinguish between departures from internal 
and external airports. For simplicity, the freeze 
horizon in Fig. 1 is represented as a circle 
centered at the arrival airport where the radius 
of the circle is the flight time from the periphery 
of the circle to the arrival airport. All airports 
inside the circle are referred to as _ internal 
airports and airports outside the circle are 
referred to as external airports. Figure 1 shows 
the freeze horizon marked by the circle Figure 1. Traffic scenario. 
centered at IAH. The two. airports — 
Austin/Bergstrom International Airport (AUS) and Dallas/Fort Worth International Airport (DFW) 
— are internal airports. The four airports- Hartsfield-Jackson Atlanta International Airport (ATL), 
Charlotte Douglas International Airport (CLT), Denver International Airport (DEN) and Miami 
International Airport (MIA) are external airports in this example. The numerical values in the 
figure indicate average flight time in minutes; for example, the average flight time from DFW to 
IAH is 44 minutes, whereas it is 135 minutes from MIA to IAH. Freeze horizon for IAH is set 
between 50 and 250 nautical-miles depending on the arrival meter-fix and the class of aircraft- 
jets, turboprops and piston-props. For jets it is set between 190 and 250 nautical-miles; this 
translates to between 46 and 60 minutes assuming an airspeed of 250 knots, which is the 
maximum allowed airspeed inside the Class-B airspace. The freeze horizon shown in Fig. 1 
corresponds to approximately 50 minutes flight time with respect to IAH. In simulations with a 
mix of CFR and non-CFR flights, only flights from DFW are assumed to be CFR flights. Flights 
from all other airports including internal airports like AUS are considered to be non-CFR flights. 

The scenario encompasses departures within 27 hours spanning 9:00 p.m. (3:00 
Coordinated Universal Time (UTC)) January 2", 2011 through midnight Central Standard Time 
(CST) (6:00 UTC) January 3", 2011 from 173 airports, including the ones in Fig. 1. Out of 808 
departures from these airports, 16 are from DFW. ACES, discussed in Section Ill, was used to 
simulate 33 hours of air traffic to enable all flights that departed within the 27 hours to land at 
IAH. Results were generated for different freeze horizon and look-ahead time values using the 
flight times computed by ACES along with the IAH arrival rate constraint in the FCFS scheduler, 
which is also described in Section III. 


Figure 2 shows the number of internal airports in this study as a function of freeze horizon in 
minutes. For the maximum value of freeze horizon of 120 minutes, 81 airports are inside the 
freeze horizon. This means that at least 92 airports (173 minus 81 airports) are outside the 
freeze horizon. All airports are external airports for a freeze horizon below 40 minutes. Six 
airports (accounting for 16 departures) are internal airports at the freeze horizon of 40 minutes. 
At the far end- 120 minutes, 417 flights depart from the 81 internal airports. DFW becomes an 
internal airport starting at a freeze horizon of 45 minutes. Note that TBFM uses a fixed freeze 
horizon in real operations. This study does not propose a changing freeze horizon; it uses 
freeze horizon as a parameter to determine the interaction between look-ahead times and 
freeze horizon values. 


B. Scheduling of CFR and Non-CFR Flights 


For CFR flights, which are only from DFW in 100 

this scenario, departure scheduling is done in 
two different ways- at gate pushback (zero look- 
ahead time) and prior to gate pushback with 
look-ahead times spanning from 15 minutes to 
60 minutes in 15-minute intervals. When 
scheduling is done at gate pushback, an arrival 
slot is assigned based on arrival time 
determined using wheels-off time and flight time 
estimates, and the assigned delay is taken on 
the ground. Taxi-out time uncertainty limited to a : 
a three-minute window of two-minutes early to 40 740 60 80 100 120 
one-minute late [-2, +1] (see Ref. 1 Section Prpezenerizpnirines) 
III.B) with respect to the nominal taxi-out time is 
permitted; because scheduling is done at gate 
pushback in this instance, the wheels-off time 
uncertainty with respect to the scheduled wheels-off time is also bounded by the [-2, +1] minute 
range of the taxi-out time uncertainty. Limited taxi uncertainty is based on the assumption that 
CFR flights are given priority on the surface during the taxi-out phase. To account for this 
uncertainty, once airborne, delayed flights can increase their nominal speed by up to 2.5% and 
early flights can reduce their nominal speed by up to 5% to get to their assigned arrival slot. This 
amount of speed adjustment was found to be adequate in this study for CFR flights to not miss 
their slots. Aircraft are thrust limited at altitude and have a very limited capability for increasing 
speed. For example, cruise Mach at cruise altitude of 35,000 feet for several Boeing 777 
variants is 0.84 and the maximum Mach is 0.89 (see Specifications by model table in Ref. 4). 
This speed difference translates to a 6% difference with respect to cruise Mach. Large deviation 
from the optimal cruise speed is also undesirable because it increases fuel consumption. 

When scheduling is done for CFR flights prior to gate pushback, an arrival slot 
corresponding to an arrival time, which is derived from estimated wheels-off time and flight time, 
is requested at a look-ahead time prior to the proposed gate departure time. This procedure 
enables earlier slot reservation. Delay assigned by the scheduler is taken on the ground. Next, 
gate departure time and taxi-time uncertainties are compared with the delay. The larger of the 
two is added to the proposed gate departure time and nominal taxi-time. The sum is determined 
to be the actual wheels-off time. Starting with a value of 1.15 for look-ahead time of 15 minutes, 
the factor for scaling gate departure time uncertainty is increased by 0.05 with increasing look- 
ahead times in steps of 15-minutes. Thus, the scale factor is 1.2 at 30-minutes and 1.3 at 60- 
minutes look-ahead times (see Ref. 5). Taxi-out time uncertainty is limited to [-2, +1] minutes. 
Once airborne, if it is determined that the delayed departures will not be able to make up for the 
lost time by increasing their speed by 2.5%, the assigned slots are freed up and made available 


Number of internal airports 


Figure 2. Number of internal airports as 
a function of freeze horizon. 
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to other flights; new slots are assigned with additional delays taken in the airborne phase by 
flying a longer route or by holding in a holding pattern. Two sets of scenarios are discussed later 
in the results section, 1) with gate departure time uncertainties scaled as a function of look- 
ahead times according to Ref. 5 and taxi-time uncertainty limited to [-2, + 1] minutes and 2) with 
controlled wheels-off time uncertainty of [-2, +1] minutes with respect to the scheduled wheels- 
off time, where the scheduled wheels-off time is obtained as the sum of proposed gate 
departure time, nominal taxi-time and scheduler assigned delay. Wheels-off time uncertainty is 
controlled by controlling the combination of the gate departure time and taxi-time uncertainties. 

Non-CFR flights from internal airports are scheduled at gate pushback based on wheels-off 
time and flight time estimates. There is no option for scheduling prior to gate pushback for non- 
CFR flights. In both types of flights, actual wheels-off time is determined based on delay 
assigned by the scheduler and taxi-time uncertainty, discussed in Section Ill. Taxi-time 
uncertainty is same for CFR and non-CFR flights; it is unique to each airport. Departures are 
allowed to recover lost or gained time, caused by taxi-time uncertainty, by increasing average 
speed by up to 2.5% or by decreasing it by up to 5% in the airborne phase. In instances when 
this adjustment is inadequate to make it to the assigned slot, the slot is freed up and made 
available to other flights; a new slot is then assigned to the flight with additional delays taken in 
the airborne phase. 

Flights from external airports are non-CFR flights that are airborne by the time they arrive at 
the freeze horizon. Upon crossing the freeze horizon, an arrival slot is assigned to them by the 
FCFS scheduler and their scheduled time of arrival at the destination airport is frozen. Any delay 
assigned by the scheduler is taken in the airborne phase. Being airborne at the time of 
scheduling means that the scheduling is unaffected by wheels-off time uncertainties, therefore 
these flights do not miss the assigned slots. 


lll. DEPARTURE DELAY MODEL AND SOFTWARE 
SIMULATION TOOLS 


Modeling of gate departure delay and taxi-out delay distributions using real-world data 
derived from the FAA’s ASPM database is described in Section III.A. Brief descriptions of ACES 
and the FCFS scheduler are provided in sections III.B and III.C. 


A. Gate Departure Delay and Taxi-out Delay Distribution Modeling 


The FAA uses ASPM data, which are accessible on the web for authorized users, for 
monitoring airport efficiency, aspects of system performance, and retrospective trend analysis 
studies. The ASPM database provides detailed data on flights to and from the 77 major U. S. 
airports (ASPM 77 airports) and flights operated by 29 major carriers (ASPM 29 carriers). 
Flights operated by ASPM carriers between ASPM airports and international and domestic non- 
ASPM airports are included. Information on airport weather, runway configuration, and arrival 
and departure rates are also available in the database. ASPM data are useful for gaining insight 
into air traffic and air carrier activity. 

To model gate departure delay and taxi-out delay distributions, hourly average gate 
departure delays and hourly average taxi-out delays in minutes were derived for each of the 77 
major U. S. airports in the ASPM database for each hour of every day in 2011 using the 
standard “Analysis-By-Airport-By-Hour Report (compared to flight plan).” Gate departure delay 
is defined as the difference between actual gate-out time and the flight-plan gate-out time. Taxi- 
out delay is defined as the difference between taxi-out time and unimpeded taxi-out time, where 
unimpeded taxi-out time is defined as the average time (a single value for each airport) taken to 
taxi out when congestion, weather, or other delay factors are insignificant. The unimpeded taxi- 
out times (nominal taxi-out times) for the airports were also obtained from the ASPM database. 


Histograms of hourly average gate departure delay and hourly average taxi-out delay in 
minutes for Hartsfield-Jackson Atlanta International airport are shown in Figs. 3a and 3b. There 
are 8,727 samples (approximately 24 hours x 365 days) in the histogram, with data missing for 
some hours on some days. Observe that data are shown from O up to 60 minutes; ASPM 
database reports non-negative delays. The maximum hourly average gate departure and taxi- 
out delays were found to be 4.6 hours and 2.4 hours, respectively. Outliers such as these are 
not shown in the figures. These histograms are modeled with log-delay distributions in this 
study. 
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To construct the log-delay distribution, zero delay values were first separated from positive 
delay values. Natural logarithm of the positive values was then computed and their distribution 
was determined. Figures 3c and 3d show the resulting histograms of the logarithm of hourly 
average gate departure and taxi-out delays, where the delays are in minutes. Both of these 
distributions are modeled by Gaussian distributions with means and standard deviation values 
computed from the log-delay distributions in spite of the fact that the log-delay distributions fail 
the Anderson-Darling test for normality. Probability Density Functions (PDF) and Cumulative 
Distribution Functions (CDF) of the log-delay and the Gaussian distributions are discussed in 
the context of modeling in the Appendix. 

To account for the zero-delay values in the model, a probability of zero-delay is computed 
as the ratio of the number of zero-delay samples to the total number of samples. For Atlanta, 
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the zero hourly average gate departure delay and taxi-out delay samples were found to be 
1,410 and 1,249, respectively. The probabilities of zero hourly average gate departure and taxi- 
out delays are thus 0.16 and 0.14 with 8,727 total samples. Data for the 77 airports are given in 
the Table 3 in the Appendix. 

The first step of simulating gate departure and taxi-out delays of a flight from one of the 77 
ASPM airports consists of generating a uniformly distributed pseudo-random number between 
zero and one. If this value is less than or equal to the zero-delay probability value in column two 
for gate departure delay in Table 3 in the Appendix and in column five for taxi-out delay in the 
same table, zero delay is assigned. When the value is greater than the zero-delay probability, 
the second step consists of generating a normally distributed pseudo-random number with 
mean and standard deviation values specified in Table 3 in the Appendix. If the numerical value 
of the pseudo-random number is x , delay 6 is obtained as 


0=e (1) 
These two steps are used to model the distributions shown in Figs. 3a and 3b. 


B. Airspace Concept Evaluation System 

ACES is a gate-to-gate simulation of air traffic developed at NASA Ames Research Center.® 
ACES simulates flight trajectories using aircraft models derived from the Base of Aircraft Data 
(BADA) (Ref. 7) and traffic data consisting of departure times and flight-plans obtained from 
recorded Aircraft Situation Display to Industry (ASDI) files. For each flight in the flight-set to be 
simulated, ACES takes as input a flight-plan, which includes names of arrival and departure 
airports, aircraft type (for example, Boeing 737), cruise speed, cruise altitude and route in terms 
of a series of latitude-longitude pairs. Typical ACES outputs include aircraft state such as 
position and velocity as a function of time, and system performance metrics such as arrival, 
departure, en-route and total delays. Validation studies in Refs. 8 and 9 have shown that ACES 
generates delays and metrics comparable to those observed in the real world. Although ACES 
can delay flights on the ground and in the air to comply with airport and airspace capacity 
constraints, this capability was not used in this study. Instead, unconstrained flight times of 808 
flights to IAH were determined without capacity constraints. These flight times and proposed 
wheels-off times are provided to the FCFS scheduler for arrival slot reservation and delay 
assignment. 


C. First-Come First-Served Scheduler 


The FCFS scheduler, as the name suggests, schedules arrival slots in the order requested 
based on airport arrival rate constraints. Input to the scheduler consists of flight times derived 
from trajectories generated during ACES simulation of unconstrained traffic and proposed 
wheels-off times obtained by summing the proposed gate departure times from flight-plans and 
unimpeded taxi-times for airports of departure from ASPM. FCFS also takes as input the taxi 
and gate uncertainty distributions described in subsection A, and the ASPM-based hourly airport 
arrival rate constraint. The airport arrival rate constraints are converted into temporal spacing 
constraints. For example, arrival rate constraints between 68 and 100 flights per hour at IAH for 
the period of this scenario translated into spacing constraints of 0.88 to 0.6 minutes between 
successive flights. Note that even when the hourly arrival capacity is greater than hourly 
demand, delays result due to competition for the same slot by multiple flights given that slots are 
spaced by at least the minimum spacing constraint described above. 

The block-diagram in Fig. 4 summarizes the procedure for scheduling and re-scheduling 
CFR flights, non-CFR flights, internal airport departures and external airport departures, which 
were discussed earlier in Section Il. Recall that the scheduler estimates the arrival time as the 


sum of proposed gate departure time, unimpeded taxi-out time and flight time; it attempts to 
assign a slot for this arrival time. If unsuccessful, it allocates the next available slot and asks for 
the flight to be delayed either on the ground or in the air depending on the state of flight at the 
time of scheduling. If a flight departs late due to uncertainties after taking the assigned delay on 
the ground, to the extent that the reserved slot cannot be reached after increasing or decreasing 
its speed, the reserved slot is reclaimed by the scheduler and a new slot and airborne delay are 
assigned. Flights from external airports are assigned an arrival slot when they cross the freeze 
horizon and the scheduling delay is taken in the air inside the freeze horizon. A new slot is 
always spaced with respect to previously assigned slots based on the airport arrival rate 
constraint, which typically varies throughout the day. From a practical standpoint, this means 
that an early slot reservation request for a slot at a later time requires an airport arrival rate 
constraint forecast at that later time. Poor arrival rate constraint forecasting can lead to the 
assignment of a slot that does not materialize at the target time or a slot that is not assigned but 
appears at the target time in the future. In this study, perfect knowledge of arrival rate constraint 
is assumed. 

The following example illustrates the scheduling procedure for a CFR flight from DFW in 
which an arrival slot at IAH is requested at 9:00 a.m. for an expected 10:00 a.m. gate departure 
time. Given that DFW unimpeded taxi-out time is 11 minutes and flight time from DFW to IAH is 
44 minutes, the scheduler attempts to reserve a 10:55 a.m. slot (arrival time estimate = 
expected gate departure at 10:00 a.m. + unimpeded taxi time of 11 minutes + nominal flight time 
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Figure 4. FCFS scheduling process. 


of 44 minutes). In this example, the 10:55 a.m. slot is already occupied therefore the next 
available slot at 10:56 a.m. and one-minute of delay are assigned by the scheduler. As a result 
of gate and taxi-out time uncertainty, the flight actually departs with a wheels-off time of 10:20 
a.m. Upon detecting wheels-off, the scheduler estimates the earliest arrival time for this flight to 
be 11:00 a.m. assuming a 2.5% average speed increase; therefore, it frees up the 10:56 a.m. 
slot at 10:20 a.m. for other potential arrivals and attempts to assign the 11:00 a.m. slot. Because 
the 11:00 a.m. slot is also unavailable, the next available slot at 11:01 a.m. is assigned. The 
flight then takes an additional minute of delay in the airborne phase to arrive at IAH at 11:01 
a.m. 


IV.RESULTS 


Results based on the procedures and tools described in sections II and Ill are presented in 
the following subsections. Subsection A describes the scenarios that were run and the 
reasoning behind them. Subsection B describes the two scenarios in which all flights are CFR 
flights as in the previous delay sensitivity study in Ref. 2. Subsection C describes the two 
additional scenarios in which both CFR and non-CFR flights are included. Lastly, Subsection D 
investigates total system delay in each scenario. These results attempt to answer the question 
of whether scheduling in advance is beneficial. All flights are IAH bound in these scenarios. In 
scenarios with CFR and non-CFR flights, only flights departing from DFW are considered to be 
CFR flights. Results are based on Monte-Carlo simulations in which each scenario is run 10,000 
times with randomized taxi time and gate departure time uncertainties. 


A. Description of Scenarios 


Three scenarios and their variations are presented in Table 1. The scenarios are 
characterized by (a) the extent of CFR flights- (1) none of the flights are CFR flights, (2) all 
flights are CFR flights and (3) only flights departing DFW are CFR flights, and by (b) wheels-off 
time prediction uncertainty- (1) wheels-off time uncertainty controlled to [-2, +1] minutes and (2) 
wheels-off time uncertainty is not controlled (sampled from distributions such as in Figures 3c 
and 3d as described in Section Ill). For the remainder of Section IV, scenarios will be referred 
by their scenario number listed in Table 1. 


Table 1. Scenario summary 


Scenario # Description Motivation 

0 None of the flights are CFR flights. Baseline scenario. 

la All flights are CFR flights and wheels-off For comparison with conclusions of 
time uncertainty is not controlled. the earlier study in Ref. 2. 

1b All flights are CFR flights and wheels-off Same reason as 1a. 
time uncertainty is controlled. 

2a Only DFW flights are CFR flights and Focal point of this study- to investigate 
wheels-off time uncertainty is not the benefits to CFR flights over non- 
controlled. CFR flights. 

2b Only DFW flights are CFR flights and Same reason as 2a. 
wheels-off time uncertainty is controlled. 


B. All Flights are CFR Flights 


This section describes the results of scenarios 1a and 1b in terms of delay per flight. These 
scenarios were run for comparison with the findings of the earlier delay sensitivity study, Ref. 2, 
in which all flights were CFR flights. 


Figure 5a presents the average weighted delay per flight in minutes, where the average is 
taken over 10,000 Monte Carlo runs and 808 departures, as a function of scheduling from 0 to 
60 minutes prior to the proposed gate departure time. Because gate departure time uncertainty 
is larger for longer look-ahead times (see Ref. 5), more flights fail to depart at the gate departure 
time assumed during slot assignment; this causes them to miss their scheduled arrival slot. 
These slots often cannot be assigned to later flights. Rescheduling to assign new arrival slots 
for missed slots also reduces the number of available slots at later times, which causes the 
overall system delay to increase. Thus, as look-ahead time increases, delay increases as 
depicted in Fig. 5a. This trend is consistent with the findings in the previous delay sensitivity 
study summarized in Ref. 2. 

Results for the second scenario- 1b are summarized in Figure 5b. Wheels-off time 
uncertainty is limited to [-2, +1] minutes in this scenario. Comparing Figs. 5b to 5a it is seen that 
limiting wheels-off time uncertainty causes the delay to remain constant with increasing look- 
ahead times and the amount of delay is the same as zero-minutes CFR look-ahead under 
uncontrolled uncertainty, which implies that there is no advantage to scheduling in advance of 
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gate departure. The main reason is that scheduling all flights equally early with small wheels-off 
time uncertainty does not change the arrival sequence and congestion at the arrival airport. 
Earlier in Ref. 2 it was determined that without gate departure time and taxi-out time 
uncertainties, scheduling in advance offered a slight advantage. Total delay could be lowered by 
re-sequencing flights and more of the scheduling delay could be passed to the ground with 
longer look-ahead times. Results in Fig. 5b do not show such an advantage because they are 
based on strict FCFS without re-sequencing. 

Although results in Figs. 5a and 5b showed no advantage to scheduling in advance of gate 
pushback, even when uncertainty is small, it is important to note that only the capacity of the 
arrival airport was taken into account by the scheduler. If other constraints such as the capacity 
of the departure airport were taken into account, it is possible that the time gained by scheduling 
in advance of gate departure could be used for improving coordination with downstream 
facilities and scheduling of surface traffic, which could lower system-wide delays. 


C. Mix of CFR and Non-CFR Flights 


Results of 2a and 2b scenarios are discussed in this section. Only flights from DFW are CFR 
flights in both these scenarios. As in Scenario 1a, gate departure time uncertainty increases for 
CFR flights with longer look-ahead times in Scenario 2a. Taxi time uncertainty is limited to -2 to 
+1 minutes with respect to the unimpeded taxi-out time. Wheels-off time uncertainty for CFR 
flights is restricted to -2 to +1 minutes in Scenario 2b like in Scenario 1b. CFR flights are 
scheduled from 0 to 60 minutes prior to their proposed gate pushback time. Non-CFR flights 
outside the freeze horizon, which is set in a range of 30 to 120 minutes of flight time to IAH, are 
scheduled when they cross the freeze horizon. These flights are not affected by gate departure 
time and taxi-out time uncertainty because they are already airborne by the time they cross the 
freeze horizon. Non-CFR flights within the freeze horizon are scheduled at their gate departure 
time; taxi-out time for these flights in Scenario 2a is obtained by sampling the taxi-out delay 
distribution and in Scenario 2b is limited to within -2 to +1 minutes. 

Results of 10,000 Monte Carlo runs for Scenario 2a are summarized in Figs. 6a and 6b. 
Figure 6a and 6b show the average weighted delay per CFR flight and per non-CFR flight, 
respectively, as a function of freeze horizon and look-ahead time. Averages were determined 
based on 16 CFR flights from DFW and 792 non-CFR flights from the other 172 airports. Figure 
6a shows that there is no advantage to scheduling in advance. The graph corresponding to 
scheduling at gate departure time is always below the graphs corresponding to scheduling at 
longer look-ahead times. Furthermore, observe the bi-modal behavior with one being the graphs 
corresponding to longer look-ahead times bunched together and the other being the graph 
corresponding to at gate departure time. Whereas the standard deviation of gate departure 
delay distribution increases by 30% from look-ahead time of zero to 60 minutes, delay does not 
change significantly in the 15-minute to 60-minute look-ahead time range; it remains primarily a 
function of freeze horizon. Figure 6b shows that the delays of non-CFR flights are very similar in 
magnitude to the delays of CFR flights obtained by scheduling at gate departure time (see zero 
look-ahead time result in Fig. 6a). This is expected because non-CFR flights inside the freeze 
horizon are also scheduled at gate departure time. Please note the different y-axis scales of 
Figs. 6a and 6b; the smaller y-axis scale in Fig. 6b is used to amplify the differences in look- 
ahead time graphs. Comparing Figs. 6a and 6b it is seen that CFR flights have a lower delay up 
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to the freeze horizon value of about 45 minutes. The reason is that CFR flights get scheduled 
prior to both external and internal non-CFR flights for freeze horizon values smaller than the 
average flight time from DFW to IAH, which is 44 minutes. For example, a CFR flight would be 
scheduled at 9:16 a.m. at the latest for an arrival slot at 10:00 a.m. With a freeze horizon of 40 
minutes or less, a non-CFR flight also competing for the same arrival slot at 10:00 a.m. will be 
scheduled at 9:20 a.m. at the earliest, after the CFR flight has already been scheduled. Delay 
increase for non-CFR flights at longer freeze horizon values seen in Fig. 6b is due to more 
airports, therefore more flights, being included inside the freeze horizon; flights from these 
airports are affected by taxi-out time uncertainty compared to those outside the freeze horizon 
that do not have uncertainty at the time of scheduling. 

Figures 7a and 7b show the Scenario 2b average weighted delay results for CFR and non- 
CFR flights, respectively. Figure 7a shows that up to a freeze horizon of 50 minutes, there is no 
advantage to scheduling early. As the freeze horizon is increased, more non-CFR flights get an 
advantage by scheduling earlier compared to CFR flights from DFW, which are about 44 
minutes away from IAH. Observe the increasing delay as a function of freeze horizon in Fig. 7a 
when CFR flights are scheduled at gate departure time compared to an initial decrease in delay 
of non-CFR flights in Fig. 7b. The graphs in Fig. 7a show the delay of CFR flights to be 
decreasing with longer look-ahead time. Since the wheels-off time uncertainty is kept small and 
independent of look-ahead time, longer look-ahead time can be thought of as a longer freeze 
horizon only for CFR flights. This enables CFR flights to schedule earlier than most non-CFR 
flights. Comparing Fig. 7a to 6a, it is seen that the delay of CFR flights is considerably lower 
when wheels-off time uncertainty is small. Smaller taxi-out uncertainty has very little impact on 
average delay of non-CFR flights; this can be observed by comparing Fig. 7b to 6b. 
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In summary, CFR flights get an advantage over non-CFR flights if they are scheduled 
earlier, which happens if the sum of their flight time to destination airport and look-ahead time is 
larger than the freeze horizon, and the wheels-off time uncertainty is small. If the gate departure 
time uncertainty is high, it is best to schedule at gate departure time. The freeze horizon 
distances of 190 to 250 nautical-miles for jets in IAH TBFM translates to freeze horizon times of 
38 to 50 minutes assuming an average speed of 300 knots. Results in Figs. 6a and 7a suggest 
that CFR flights from DFW to IAH should be scheduled at gate departure time for this range of 
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freeze horizon times; there is no advantage to scheduling earlier even with controlled 
uncertainty. 


D. Total System Delay 


Delays for CFR and non-CFR flights were examined in subsections B and C for the five 
scenarios listed in Table 1; in this section, average system delay measured as the sum of arrival 
delays (includes both ground and airborne delays) of all flights to IAH is examined. 
Combinations of look-ahead times and freeze horizons were examined to determine maximum 
and minimum average total delays. The averages are again based on 10,000 Monte Carlo runs. 
Results are summarized in Table 2 and in Fig. 8. Minimum and maximum average total delays 
were found to be 13.8 hours for a freeze horizon of 30 minutes and 31.2 hours for a freeze 
horizon of 120 minutes, respectively, for the baseline scenario without CFR flights (Scenario 0 in 
Table 1). Average total delays for the baseline scenario for the range of freeze horizon values 
are shown in Fig. 9a in the Appendix. The least amount of delay is achieved in this scenario, 
compared to the delays in other scenarios listed in Table 2, because all flights are scheduled in 
the order of their arrival. One of the findings in the previous study in Ref. 2 is that system delay 
increases as arrival slots for more flights are scheduled out of the nominal arrival order. 


Table 2. Total system delay summary 


Scenario # | Delay (Hrs.) | Look-ahead time and freeze horizon parameters 
0 (min) 13.8 Look-ahead time: all, Freeze horizon: 30 
0 (max) 31.2 Look-ahead time: all, Freeze horizon: 120 
1a (min) 40.2 Look-ahead time: 0, Freeze horizon: all 
1a (max) 90.0 Look-ahead time: 60, Freeze horizon: all 
1b 40.2 Look-ahead time: all, Freeze horizon: all 
2a (min) 14.2 Look-ahead time: 0, Freeze horizon: 30 
2a (max) 34.0 Look-ahead time: 60, Freeze horizon: 120 
2b (min) 14.2 Look-ahead time: 0, Freeze horizon: 30 
2b (max) 30.4 Look-ahead time: 60, Freeze horizon: 120 


Because there are no internal airports within the 30-minute freeze horizon, arrival slots are 
scheduled in the order they would have naturally arrived at IAH. Additionally, because all flights 
are airborne by the time they are scheduled at the freeze horizon, their scheduled arrival slots 
are unaffected by gate departure and taxi- 
out delays, which cause flights to be 
rescheduled in the other scenarios. 

The maximum average system delay 
occurs in Scenario 1a with a look-ahead 
time of 60 minutes when all flights are CFR 
flights and gate departure time uncertainty 
is not controlled. Results are independent saul jal 
of freeze horizon as seen in Fig. 10a. This 56 i) '. ' ' 
is to be expected because there is no 46 1 
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scheduled slots and (2) all flights are scheduled in the order of departure, i.e. when they leave 
the gate. When flights miss their arrival slot, they have to be rescheduled and additional delay 
has to be taken while airborne. Scheduling in the order of departure leads to higher delay 
compared to scheduling in the order of arrival. 

The next set of results for Scenarios 2a and 2b in Table 2 are similar to the results for 
Scenario 0 as expected; Scenario 2 has only 16 fewer non-CFR flights (16 CFR flights from 
DFW and 792 non-CFR flights from other airports) compared to Scenario 0 with 808 non-CFR 
flights. Graphs in Fig. 12a in the Appendix show the average delay to be mostly a function of 
freeze horizon. Delay increases a bit when scheduling is done prior to gate departure. Graphs 
for Scenario 2b in Fig. 13a in the Appendix show the average delay to be almost independent of 
look-ahead time. This is expected because wheels-off time uncertainty is limited. 

In summary, the minimum values for scenarios 0, 2a, and 2b in Table 2 show that 
scheduling in the order of arrival, which is accomplished by setting the freeze horizon to 30 
minutes, leads to small average system delay. Maximum values for the same scenarios in Table 
2 show that longer freeze horizon leads to large average system delay. Most system delay 
occurs when all flights are CFR flights and departure time uncertainty is not controlled. 


V. CONCLUSIONS 


An earlier delay sensitivity concluded scheduling should be done at gate pushback because 
of gate departure time and taxi time uncertainties. It also concluded that there is a slight 
advantage of scheduling a flight early for an arrival slot in that the delay needed for temporal 
spacing can be taken on the ground as opposed to in the air, provided the wheels-off time 
uncertainty can be controlled to a small value (like -2 minutes early to +1 minute late). All flights 
were Call For Release flights in the earlier study. 

The current study is based on the same simulated air traffic data as the earlier one. Five 
different scenarios with 808 Call For Release and non-Call For Release flights from 173 airports 
to Houston George Bush Intercontinental Airport (IAH) were examined. These scenarios either 
assumed gate departure time uncertainty that increased with look-ahead time or wheels-off time 
uncertainty that could be limited to -2 minutes early to +1 minute late with respect to the 
scheduled wheels-off time. Average weighted delays (sum of one-time ground delay + two- 
times airborne delay) and average arrival delays were computed for these scenarios based on 
10,000 Monte Carlo runs using the First-Come First-Served method for scheduling arrivals. 

Results show that increasing gate departure time uncertainty as a function of look-ahead 
time causes the weighted delay to increase as a function of look-ahead time when all flights in 
the scenario are Call For Release flights. It is detrimental to schedule early in this scenario. It 
also follows that there is no advantage to early scheduling even if wheels-off time uncertainty is 
small in this scenario because every Call For Release flight can schedule equally early. These 
findings are in agreement with the earlier study. 

Analysis of scenarios with 16 Call For Release flights from Dallas-Fort Worth Airport (DFW) 
and 792 non-Call For Release flights from 172 other airports showed that Call For Release 
flights do get an advantage over the non-Call For Release flights from airports inside the freeze 
horizon if the sum of the flight time to the destination airport and look-ahead time is longer than 
the freeze horizon and the wheels-off time uncertainty is small. If gate departure time 
uncertainty is large, Call For Release flights should be scheduled at gate pushback. 

Results in this report show that Call For Release flights from Dallas-Fort Worth to Houston 
George Bush Intercontinental do not benefit by scheduling prior to gate departure time for a 
typical Traffic Management Advisor freeze horizon setting of 38 to 50 minutes or an equivalent 
distance of 190 to 250 nautical-miles. From a system delay (sum of arrival delays of all flights) 
perspective, the order of scheduling was found to be important. Results show that scheduling in 
the order of arrival leads to less delay compared to scheduling in the order of departure. 
Increasing the freeze horizon leads more flights to be scheduled in the order of departure, which 
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increases system delay. The largest amount of system delay was accrued at the longest look- 
ahead time in the scenario with all Call For Release flights primarily due to increasing gate 
departure time uncertainty with longer look-ahead time. 


Table 3 


APPENDIX 


shows the parameters for modeling gate departure and taxi-out time delay 


distributions shown in Figs. 3a and 3b for the 77 ASPM airports following the procedure 
described in Section IIl.A. 
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Table 3. Hourly average gate departure and taxi-out delay model data 


Hourly Average Gate Departure Delay Hourly Average Taxi-out Delay 


Airport 
Code 


Zero- Weak Standard Zero- Klean Standard 
Delay Loo-bela Deviation Delay LeoDela Deviation 
Probability | -°9°-'©'*Y | Log-Delay _| Probabilit gveley | Log-Dela 


ATL 
AUS 
BDL 
BHM 
BNA 
BOS 
BU 
BUR 
BWI 
CL 
CL 


DA 
DA 


DE 
DFW 
DTW 
EWR 
FL 


1.089 
0.966 
0.842 
1.280 
0.934 
0.870 
1.099 
0.772 
1.264 
1.036 
1.222 
0.907 
0.928 
0.666 
1.144 
0.901 
0.914 


0.938 0.186 


HNL 


0.843 


0.939 


HPN 


0.966 1.300 


IND | 0.204 | 1.846 | 1163 | 0.249 | 0057 | 1216 


JAX 


Table 3. Hourly average gate departure and taxi-out delay model data (contd.) 


Hourly Average Gate Departure Delay Hourly Average Taxi-out Delay 


Airport 
Code 


JFK 
LAS 
LAX 
LGA 
LGB 
MCI 
MCO 
MDW 
MEM 
MHT 
MI 
MK 
MS 
MSY 
OA 
OGG 
OMA 


ON 


OXR 


PD 
PH 
PH 
PIT 
PS 
PV 
RDU 
RFD 


Zero- Mean Standard Zero- Maan Standard 
Delay loadela Deviation Delay eqns Deviation 
Probability | -°9°-°'?¥ | Log-Delay _| Probabilit gve'ay | Log-Dela 


1.048 
0.732 
0.883 
1.276 
1.286 
0.699 
1.123 
0.988 
0.982 
1.061 
0.964 
1.207 
1.184 
0.797 
0.459 
1.000 
0.863 
1.147 
0.826 
1.215 
1.074 
1.206 
0.926 
0.595 


1.074 


SAN 
SAT 
SDF 
SEA 


0.792 
0.909 


SJC 
SJU 


Table 3. Hourly average gate departure and taxi-out delay model data (contd.) 


Hourly Average Gate Departure Delay Hourly Average Taxi-out Delay 


Airport | 5 .¥0-Delay | Mean Standard Zero- Mean Standard 
Code | probability | Log-Delay | Deviation | Delay | | ooh atay | Deviation 
Y pog'8¥ | Log-Delay _| Probability | -°9°-©"®Y | Log-Delay 


SLC 
SMF 


STL 
TEB 
TPA 
TUS 
VNY 


Figures 9 through 13 present ground, airborne and weighted delay results for CFR and non- 
CFR flights, and system delay results for the five scenarios discussed in Table 1. Each figure is 
composed of sub-figures ‘a’ through ‘m’. The contents of the sub-figures are listed in Table 4. 
Color coded graphs in the sub-figures are for look-ahead times from at pushback (zero minutes) 
to 60 minutes prior to expected gate pushback time. The average delay values are shown as a 
function of freeze horizon in the graphs. Figure 9 shows the results for the baseline scenario 
without CFR flights (Scenario 0 in Table 1). Results for Scenario 1a and 1b (see Table 1) with 
all CFR flights are given in Figs. 10 and 11. Mixed scenarios- 2a and 2b results with CFR flights 
from DFW are provided in Figs. 12 and 13. 


Table 4. Contents of sub-figures 
Sub-figure Content 
Average system delay (total arrival delay at IAH) 
Average weighted delay per CFR flight 
Average ground delay per CFR flight 
Average air delay per CFR flight 
Average weighted delay per non-CFR flight 
Average ground delay per non-CFR flight 
Average air delay per non-CFR flight 
Average weighted delay per non-CFR external flight 
Average ground delay per non-CFR external flight 
Average air delay per non-CFR external flight 
Average weighted delay per non-CFR internal flight 
Average ground delay per non-CFR internal flight 
Average air delay per non-CFR internal flight 
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Figure 9. Baseline scenario without CFR flights (Scenario 0 in Table 1). 


system delay 
weighted delay/CFR 


air delay/non-CFR I. ground delay/non-CFR internal 
weighted delay/non-CFR external m._ air delay/non-CFR internal 
ground delay/CFR ground delay/non-CFR external 

air delay/CFR air delay/non-CFR external 

weighted delay/non-CFR ~_k. weighted delay/non-CFR internal 

ground delay/non-CFR 


SG 


2 < < 
4 
20 =8 = = 
“15 g oa Bs 
+ 30 Qo o Dm 
« 45 a) Ow a 
60 5 E 5 co 
=e S o> ooganopnenpaneqaean! BE 
Legend: CFR Look 40 80 120 40 80 120 
Ahead (min) Freeze Horizon Freeze Horizon Freeze Horizon 
a. ; Cc. 
= eo =) x 
eis Hyg XO OOK EK DCM ERIK IK = a Pe 
wo = 2 =o =o 
3 ee tt a o 3 co a co 
Cc Cc = Cc 
oe ool oo i 
=o =i T =" = * 
40 80 120 40 80 120 40 80 120 40 80 120 
Freeze Horizon Freeze Horizon Freeze Hotzon Freeze Horizon 
F e. i g. 
eo <o eo - eo _- 
EE ET —E° —E~ 
= Bo - Be 
=O == = = 
ac 8s 8s as 
¢ 4 4 i= 
gc gc gc gc 
=* aot T =* T =* T 
40 80 120 40 80 120 40 80 120 40 80 120 
Freeze Horon Freeze Horizon Freeze Horizon rraeze Horzon 
B i. j. . 
£°e €°o 
Ew Ew 
s s 
=O =o 
as ac 
& So. 
=v =* 
40 80 120 40 80 120 
cia ial Freeze Horizon 
, m. 


Figure 10. Results for Scenario 1a in which all flights are CFR flights and wheels-off time uncertainty 
is not controlled. 
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Figure 11. Results for Scenario 1b in which all flights are CFR flights and wheels-off time uncertainty 
is controlled. 
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Figure 12. Results for Scenario 2a in which only DFW flights are CFR flights and wheels-off time 
uncertainty is not controlled. 
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Figure 13. Results for Scenario 2b in which only DFW flights are CFR flights and wheels-off time 
uncertainty is controlled. 
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Modeling of Gate Departure Delay and Taxi-out Delay Distributions with Normal 
Distributions 


Earlier in Section IIl Subsection A, it was stated that the log of the gate departure delay 
and log of the taxi-out delay distributions such as the ones in Figs. 3c and 3d for the 77 ASPM 
airports, listed in Table 3 in this Appendix, were modeled as Gaussian distributions with mean 
and standard deviation values given in Table 3. It turns out that the log-delay distributions, while 
they appear to be symmetric, are not Gaussian distributions based on the Anderson-Darling test 
for normality. Given this fact, is it reasonable to model the log-delay distributions with the 
Gaussian distributions? This assessment was made by comparing the Probability Density 
Functions (PDF) and Cumulative Distribution Functions (CDF) of the log-delay distributions 
derived from the histograms such as the ones in Figs. 3c and 3d with those derived from 
analytical formulas for Gaussian distributions. 

Figure 14 and 15 show the PDF and CDF of the ATL gate departure log-delay distribution 
and the Gaussian distribution with the same mean and standard deviation as that of the log- 
delay distribution. PDF for the log-delay distributions is computed from the discrete histogram as 
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Figure 14. Comparison of ATL gate 
departure log-delay and 
Gaussian PDFs. 


Figure 15. Comparison of ATL gate 
departure log-delay and 
Gaussian CDFs. 


follows. Let n, be the number of samples in bin i of size Ax, of the log-delay histogram with / 
number of discrete bins. The probability density is then 


i= = (2) 
lsi<I 
where /, is the probability density associated with bin 7. It follows from Eq. (2) that 
yeas 
2. f,Ax, = 7 =1 (3) 
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The probability p, of log-delay value being less than or equal to the value represented by bin 7 
is computed by summing the probability density up to bin 7. Thus, 


R=) Ae (4) 


l<k<i 
The analytical formula for probability density for continuous Gaussian distribution is 


! oe 
f(x) — (=e (5) 


where yw is the mean and o is the standard deviation. The probability of being less than or 
equal to x is obtained by integrating the probability density defined in Eq. (5) as follows 


p(x) = | f(v)dy (6 
It can be shown that 
p(x) = Ak ve =F) a 


where erf(x) is the error function defined as 


ae ee 
af(e)=——[e dy (8) 
0 


Figure 14 shows the probability densities 
computed using Eqs. (2) and (5), respectively, 
and Fig. 15 shows the probabilities computed 
using Eqs. (4) and (7), respectively. The 
difference of log-delay probability with respect to 
the Gaussian probability is shown in Fig. 16. 

To improve the Gaussian distribution model of 
the log-delay distribution, either the standard 
deviation or both the mean and the standard 
deviation of the Gaussian distribution can be 
determined by either linearizing the PDF or the 
CDF of the Gaussian distribution and minimizing 


Probability difference 


the square of the difference with respect to the -4 9 7 : 6 
log-delay distribution PDF or CDF. Log hourly average gate departure delay (minutes) 
Method I: Linearizing the PDF Figure 16. Absolute difference of ATL 


Let, 4 and o be the mean and the standard gate departure log-delay CDF 
ne = : with respect to Gaussian CDF. 
deviation values obtained from the log-delay 


distributions. Nominal probability density is obtained using Eq. (5) as 


I (x)= [—- > (9) 


Residual (error) with respect to the probability density given by Eq. (2) is 


df=f-f (10) 
Rewriting Eq. (5) 
2 

eps) eco (1) 
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Linearizing with respect to mean and standard deviation 
we du (z |= df (12) 

Ou O 


Differentiating Eq. (11), evaluating the partial derivatives about nominal mean and standard 
deviation values, and collecting terms, Eq. (12) can be written in the matrix notation as 


(x~ w) (x= mW) - 0" § [du 
2 7 A 3 =f (13) 
— 26 = oO do 
Equation (13) is of the form 
AX =B (14) 
Least-squares solution for the unknowns is given as 
Pa 
BY) =X=(A'AA'B (15) 
do | 


where 4’ is the transpose of matrix 4 and (4'A)"' is the inverse of (4'A). The mean and 
standard deviation of the Gaussian distribution model can be determined as 


(16) 


M=u+du 
o=o+do 

Figures 17 and 18 show the comparison of the ATL gate departure log-delay and Gaussian 
PDFs and CDFs, where the Gaussian PDF and CDF are obtained with the mean and standard 
deviation values obtained using Eq. (16). The absolute difference of the probability of the log- 
delay distribution with respect to the Gaussian distribution is shown in Fig. 19. Comparing Fig. 
19 to Fig. 16 it is seen that the maximum difference between the log-delay CDF and the 
Gaussian model CDF based on the least-squares estimates in Eq. (16) is reduced from 7% 
(with respect to the maximum possible difference of 1 between the two CDFs) to 3% in this 
example. This however was not found to be true in other examples. Keeping the mean of the 
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Gaussian same as the mean of the log-delay distribution and a least-squares estimate of the 
standard deviation (using the second term in A matrix in Eq. (14)) was found to be robust in 
reducing the CDF difference. 
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Gaussian PDFs (Method 1). Gaussian CDFs (Method 1). 


Since the objective is to determine the 
mean and the standard deviation values that 
minimize the difference between the CDFs of the 
log-delay distributions and the Gaussian 0.06 Lash less hasthian aise boaantondavandsatacasoca satin doideiea deo sae annke an deATE RRS 


distribution models, an alternative least-squares 8 g 
estimation can be accomplished based on = & 

linearization of the CDF in Eq. (7). — O. 

2 0. 

Method II: Linearizing the CDF 8 : 

Let, y be a transformed variable defined as a 


my 


y= 
oV2 


Equation (7) can be rewritten in terms of y as 
follows 


-4 z 
Log hourly average gate departure delay (minutes) 
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Figure 19. Absolute difference of ATL 
gate departure log-delay 
CDF with respect to 
Gaussian CDF (Method 1). 


1 
Py) = sll +erfO)] (18) 
Partial derivative of p with respect to y is obtained as 


Op _ 1 Oerf(y)) _ 1 oe 


= (19) 
dy 2 oy Vx 


Linearizing Eq. (18) with respect to the mean and the standard deviation 


where 
dp=p-p 
Using Eqs. (19) and (17) 
X—-H 7 
Op _ Op oy _ ws) 
Ou dyOu oV2n 
and 
ap _ ey («-w) (74) 


These two equations can be written in the form of Eq. (14) with 


1) 
A= : eee 
ov2n 
B=dp 
d 
ole 

do 


The least-squares solution is then obtained via 
Eq. (15) with the matrix and vectors defined in 
Eq. (24). The mean and the standard deviation 
values for the Gaussian model are then 
obtained using Eq. (16). 

Figures 20 and 21 show comparisons of 
the ATL gate departure log-delay and Gaussian 
model PDFs and their CDFs, respectively, 
where the Gaussian model is based on the 
mean and standard deviation values obtained 
via Eqs. (24) and (16). Figure 22 shows the 
difference between the two CDFs in Fig. 21. 
Comparing Fig. 22 to Fig. 19, it is seen that the 
maximum of the absolute difference between 
the log-delay CDF and the Gaussian CDF is 
reduced further by 1% from 3% to 2%. The fit 
with respect to the PDF is a bit worse as 
observed comparing Fig. 20 to Fig. 17. The 
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Figure 20. Comparison of ATL gate 
departure log-delay and 
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(21) 


(22) 
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Gaussian PDFs (Method II). 


least-squares procedure based on CDF is the preferred procedure for both gate departure delay 
and taxi-out delay distributions. The revised mean and standard deviation values for the 
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Gaussian models for gate departure delay and taxi-out delay for the 77 ASPM airports are listed 
in Table 5. 
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Figure 21. Comparison of ATL gate Figure 22. Absolute difference of CDFs 
departure log-delay and in Fig. 21 (Method Il). 


Gaussian CDFs (Method Il). 


While the conclusions of this report are not expected to change because of the revised 
mean and standard deviation values, the Gaussian distributions constructed using the values in 
Table 5 are closer representations of the log-delay distributions derived from the real air traffic 
data. In addition to the mean and standard deviation values determined using Eqs. (24) and 
(16), Table 5 lists the maximum absolute difference between the log-delay CDFs and the 
modeled Gaussian CDFs. Observe that the maximum absolute difference for gate departure 
delay CDFs is 6% for the 77 airports listed in Table 5. This difference for taxi-out delay is less 
than 10% except for Gary/Chicago International (GYY), Oxnard (OXR), Chicago Rockford 
International (RFD), Stewart International (SWF), Teterboro (TEB) and Van Nuys (VNY); values 
for these six airports are marked in red. The number of samples for constructing Gaussians for 
these airports were small in comparison to the 8,727 samples for Atlanta, mentioned earlier in 
Section Ill Subsection A. GYY, OXR, RFD, SWF, TEB and VNY had 65, 176, 306, 1245, 1239 
and 498 samples, respectively, in the 2011 data derived from FAA’s ASPM database employed 
in this report. These samples should be augmented with samples from other years from the 
ASPM database to construct reliable Gaussians for modeling taxi-out log-delay distributions for 
these airports. 
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Table 5. Hourly average gate departure and taxi-out log-delay Gaussian models 


Airport 
Code 


ATL 
AUS 
BD 
BHM 
BNA 
BOS 
BU 
BUR 
BWI 
fen 
GE 
CVG 
DA 
DAY 
DCA 
DE 
DFW 
DTW 
EWR 
FL 
GYY 
HN 
HOU 
HP 
IND 
JAX 
JFK 
LAS 
LAX 
LGA 
LGB 


OUTY hia ee Benes Hourly Average Taxi-out Delay 


Maximum Maximum 
seal Absolute CDF ee Absolute CDF 
error (%) error (% 


1407] 1186] | 517} 898} 
1430] 1.015) | 0.082] 1-416 | 
_ 1670] 1.096} 4 | 0.204] 0.920 | 
_ 1682] 1.059) 4 | 0.438 | 0.858 | 
_ 1938] 1.208) 4 | 89] 1.280 | 
1506] 1.286) 4 | 0.304] 0.859 | 
3.034] 0.951] 4 | -0.384 | _0.008 | 
_ 1.078] 0993] 4 | 0.268] 0.814 | 
_ 1943] 0953] 4 | 0.098 | _ 0.898 | 
2611] 0847] 4 | 0.223] 1.292 | 
_ 1731] 0740] 4.971] 0.615 | 
1887] 0955] 21.922) 0.947] 


1630] 1336f 3] 0454] 136] 


Table 5. Hourly average gate departure and taxi-out log-delay Gaussian models (contd.) 


HOuny he ea Departure Hourly Average Taxi-out Delay 


Airport Maximum Maximum 
Cod 
eee Se ee Absolute CDF haas ie Absolute CDF 
error (%) error (% 


Mcl__|_1633| 1.058} 3] 0281) 49 
Mco_| 1.633] 0.963] 3 0.961] 0.622] | 
Mow | 2.052] 0849/3] | 808} 
MEM_|_2.115| 1.045) 3] 0.330) 1.120) 
MH 
MIA 
MK 
MS 
MSY | 1463] 1.245/ 4] 0.316 | 0.924 
OAK | 1.761} 4.001) 4] 0.026} 1.052 
oGG_| 0998] 1188] 4] 0.002] 02.951 
OMA_| 1.699} 4.153] 4] 0.147] 1.150 
ON 
OxR | 2.989] 0.999] 4] 1.153] 0.230 
2.190] 1052/4] 0.328] 0.953 
PDX 
PH 
PHX 
PIT 
PSP 
PvD__| 1858] 1.239) 4 0.3897 | 1.174 
RD 
RFD | 2414/ 1315| 6 | 0.550] 0.240 
Rsw_ | 1512/1258) 4 | 0.325 | 0.974 
SAN 
SAT | 1.573] 4.108] 4] 0.582] 0.750 
SDF__| 1.953] 1187/4] 0.036] 1.239 
SEA 
SFO 
SJC_ | 1453] 4.131] 4 91] 0.852 
SJU 
SLC 
SMF 
SNA__| 1496] 1106] 4] 67] 1028} 
stu | 1.779} 0.946] | 112] 0-850 


Table 5. Hourly average gate departure and taxi-out log-delay Gaussian models (contd.) 


EOutly fake Pepantte Hourly Average Taxi-out Delay 


Apert Maximum Maximum 
Cod 
oce | Mean Slate Absolute CDF | Mean Saath Absolute CDF 
error (%) error (% 


swe | 2631 1.285{ |S 1629 10 
TEB | 3.294] 0671] S| 0.864] 0.182] 1 | 
TPA 
TUS 
VNY__| 3.280] 0.790] 4] 0.545] 0.045] 14 | 
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